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PACS 47 . 27 . Jv - High-Reynolds- number turbulence 

PACS 47 . 27 . Gs - Isotropic turbulence; homogeneous turbulence 

PACS 47 . 80 . Cb - Velocity measurements 

Abstract. - We report experimental results on the relative motion of pairs of solid spheric 
particles with initial separations in the inertial range of fully developed turbulence in water. The 
particle densities were in the range of 1 J$ p P /pf ~ 8, i.e., from neutrally buoyant to highly 
inertial; and their sizes were of the Kolmogorov scale. For all particles, we observed a Batchelor 
like regime, in which particles separated ballistically. Similar to the Batchelor regime for tracers, 
this regime was observed in the early stages of the relative separation for times t ;S O.lto with to 
determined by the turbulence energy dissipation rate and the initial separation between particle 
pairs. In this time interval heavier particles separated faster than fluid tracers. The second order 
Eulerian velocity structure functions was found to increase with density. In other words, both 
observations show that the relative velocity between inertial particles was larger than that between 
tracers. Based on the widely used, simplified equation of motion for inertial point-particles, we 
derived a model that shows an increase in relative velocity between inertial particles. In its scale 
dependence, however, it disagrees quantitatively with the experimental results. This we attribute 
to the preferential sampling of the flow field by inertial particles, which is not captured by the 
model. 



Many natural and industrial phenomena involve the in- 
teraction between turbulent flows and inertial particles, 
i.e., particles that do not passively follow the fluid mo- 
tion. The transport of inertial particles by turbulent flows 
plays a fundamental role in a vast range of systems, such 
as, sedimentation in estuaries and rivers [JJ, the dynamics 
of plankton in the ocean [HE], deep sea land slides [I], 
dust in tornadoes [5], sandstorms over deserts or on Mars 
[5] , the dynamics and collisions of water droplets in clouds 
[7] and to the dynamics and clustering of interstellar dust 
in planet formation [5] . Quantitative experimental studies 
on the dynamics of inertial particles in turbulence there- 
fore provide much needed data to verify theoretical models 
and to test both simulations of idealized particle dynamics 
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in turbulent flows [SHU], as well as, coarse grained simu- 
lations of natural and technical flows that rely on so called 
"subgrid" models that parametrize the small, unresolved 
scales of turbulence [H] . 

The dynamics of single inertial particles in turbulence 
has been studied experimentally since the pioneering work 
by Snyder & Lumley [T3] and has been a topic under exten- 
sive experimental [TH - H7] . theoretical [THHH] and numeri- 
cal [MTU] investigations. Interested readers are referred to 
a recent review |20j . 

In this letter, we report an experimental investigation of 
relative dispersion between two inertial particles with ini- 
tial separations in the inertial range of fully developed tur- 
bulent water flow. We focused on solid spherical particles 
with densities larger than or nearly the same as the fluid. 
The size of particles were that of the Kolmogorov scale, 
the smallest scale in turbulence. The particle Stokes num- 
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Fig. 1: (a)The experimental setup, (b) Schematic side view of 
the apparatus. The solid lines represent the walls of the glass 
container, the dashed lines represent the removable cover. Dur- 
ing the experiments, the "aquarium" was filled with water to 
a level above the cover, as shown by the grey shadow. Two 
propellers created a highly turbulent flow in the hexagonal- 
cylindrical shaped region determined by the cover and the 
aquarium walls, while the water above the cover remained qui- 
escent. 



ber (see definition later) of the particles was in the range 
between 0.09 and 0.5. For these particles, we observed a 
Batchclor like regime, in which particles separated ballis- 
tically. Similar to the case of fluid tracers this 
regime existed in the early stage of relative separation for 
times below a timescale to determined by the turbulence 
energy dissipation rate and the initial separation between 
particle pairs. In this Bachelor like regime, particles with 
larger Stokes numbers separated faster, i.e., the relative 
velocity between heavy particles were larger than that be- 
tween fluid tracers at the same separations. This was also 
supported by the measured second order Eulcrian velocity 
structure functions. Based on the widely used, simpli- 
fied equation of motion for inertial point particles |18j . 
we derived a model that captures the observed increase 
in relative velocity between inertial particles. In its scale 
dependence, however, it disagrees quantitatively with the 
experimental results. This can most likely be attributed 
to the preferential sampling of the flow field by inertial 
particles, which is not captured by the model. 

We studied the motion of solid spherical particles in a 
von Karman swirling water flow. The turbulent flow was 
generated by two counter-rotating propellers submerged 
in an aquarium and was confined on top by a removable 
cover with rubber seals along the rim (see Figure [T]). Two 
underwater air-motors drove the baffled propellers so that 
a highly intense turbulent flow could be produced in an 
apparatus with moderate size. The diameter of the pro- 
pellers was 28 cm. The turbulence chamber, shaped as a 
hexagonal cylinder, measured 40cm along the axis of the 
propellers and 38cm in both height (vertically) and width 
(horizontally) in the cross-section. The rotating axis of the 
propellers was in the horizontal direction so that the heavy 
particles settling down towards the bottom of the appa- 
ratus were entrained by the strong sweeping of the fluid 
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75 ±8 


1.04 


0.27 
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7.8 


75 ± 15 


1.04 


0.50 



Table 1: Characteristics of the different particles used in the ex- 
periments. The Stokes number is defined as St = jjp(d p /r]) 2 , 
where d p is the particle diameter, /3 = 3p//(2p p + p/) is the 
modified density ratio with p p and p/ being, respectively, the 
particle and the fluid densities, and r\ is the Kolmogorov length 
scale of the turbulence. 

near the bottom surface. In this way particles stayed sus- 
pended in the measurement volume. Measurements were 
conducted at the center of the apparatus where the mean 
fluid velocity was small [53] . 

We used three types of particles with the same average 
size, but different density, as shown in Table [TJ We ad- 
justed the motors speed such that the Kolmogorov scale 
of the flow was close to the particle size. In this case, the 
almost neutrally buoyant polystyrene particles, behaved 
as tracer particles, as shown in previous experiments [25] . 
The other two types of particles: the Barium-Titanium 
glass spheres and the stainless steel spheres, are signifi- 
cantly heavier than water. To characterize particle inertia, 
we define the Stokes number as St = jjp(d p /ri) 2 , where 
dp is the particle diameter, (3 = 3pf/(2p p + pf) is the mod- 
ified density ratio with p p and p/ being, respectively, the 
particle and the fluid densities, and r\ is the Kolmogorov 
length scale of the turbulence. The Stokes numbers for 
these particles were 0.09, 0.27 and 0.50, respectively. 

We measured three-dimensional particle trajectories 
with high spatial and temporal resolutions using La- 
grangian Particle Tracking [55][53] with three high speed 
CMOS cameras (Phantom V10, manufactured by Vi- 
sion Research Inc., Wayne, USA). The particle velocities 
and accelerations of the particles were then obtained by 
smoothing and differentiating the trajectories [35]. The 
measurement volume of LPT was determined by the in- 
tersection of the fields of view from all three cameras. To 
avoid the effect of biased sampling due to the shape of the 
measurement volume, we used only particle trajectories in- 
side a sphere with diameter of 22 mm located at the center 
of the apparatus. We measured 3 x 10 7 data points for the 
polystyrene particles, 8 x 10 7 for the glass particles and 
3 x 10 6 for the steel particles. The LPT provided simul- 
taneous measurements of multiple particles and hence the 
data could also be used to extract Eulerian statistics even 
though the LPT technique itself is inherently Lagrangian. 
For example, knowing velocities of different particles at the 
same instant, i.e., on the same "frame", one could mea- 
sure velocity increments at distances determined by par- 
ticle positions. The Eulerian velocity structure functions 
were then obtained by collecting the statistics of velocity 
increments over many frames. The turbulence properties, 
such as the energy dissipation rate per unit mass e and the 
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Fig. 2: Illustration of the different methods used to measure 
the energy dissipation rate e (see text). The horizontal, bold 
dashed-line shows the measured dissipation e = 0.037 m 2 /s 3 , 
the thin dashed-lines represent ±5% of this value. The vertical 
line corresponds to the separation r = 407?. 



integral scale L, were inferred from the Eulerian statistics 
using polystyrene particles as fluid tracers. In the iner- 
tial range 77 <C r <C L of homogeneous and isotropic tur- 
bulence, neglecting intermittency corrections, the second 
order longitudinal and transverse velocity structure func- 
tions should scale as 



and 



D LL (r) =C 2 {erf/\ 



Dnn(t) = gCMerf/ 3 , 



(1) 



(2) 



where C2 is expected to be a universal constant. We used 
C2 = 2.1 as suggested from a compilation of available 
data [26] . In addition, there were two exact inertial range 
relations: the Kolmogorov's celebrated "four-fifth law" : 
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and a theorem on the velocity-acceleration mixed struc- 
ture function I27H301: 



(<5 r u ■ S r a) 



-2e. 



(4) 



We measured e using each of the four equations above. 
For length scales in the inertial range (r/77 2; 50) we found 
that measurements from all four methods were in agree- 
ment within ±5%, as shown in Figure [U We took the 
weighted average from all four measurements as e. The 
integral length scale was estimated as L = u' 3 /e, where 
u' was the fluctuating velocity. Other parameters, such 
as the Kolmogorov length and time scales, were then ob- 
tained from their standard definitions. The turbulence 
properties and other experimental parameters are sum- 
marized in Table [5] 
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Table 2: Parameters of the experiment, u is the root-mean- 
square of the velocity, e is the energy dissipation rate per unit 
mass. L = u' A /e is the integral length scale. 77 = (^ 3 /e) 1//4 
and tv, = (i//e) 1//2 are the Kolmogorov length and time scales, 
respectively, where v is the kinematic viscosity of the fluid. Nf 
is the frame rate of the camera, in frames per r, ? , and — is 
the resolution of the camera compared to 77. 



We now consider the separation of two particles in time 
R(i) = ra(t) — ri(i) (where i"i(i) stands for the position 
of particle i time t). For fluid tracers, it is well known 
(see e.g. [31]) that particle velocities are uncorrelated at 
larger distances (R ^> L) and therefore particles separate 
diffusively (R oc i), while they separate chaotically at sepa- 
rations below Kolmogorov scale (R <C rj). The interesting 
question is on the relative dispersion of particles with sep- 
arations in the inertial range, i.e., 77 <C R L, or equiv- 
alently 77 < R < L and t <C T L , where R a = |R(t = 0) 
is the initial separation and Tj, = L/u' = (i 2 /e) 1 ^ 3 is 
the large eddy turnover time. Theoretical work [32,33 
suggests the following two regimes for tracer dispersion: 



(SR ■ SR) = 



(6MR0) ■ 5 r v(R ))t 2 = 

for t < to = (i?oA) 1/3 

and 

get 3 for t <C t < T L 



(5) 

where SR(t) = R(t) — R(t = 0) is the vectorial separation 
increment, to may be regarded as the life time of an eddy 
of scale Rq, and g is a dimensionless coefficient, known 
as the Richardson constant, expected to be universal and 
independant of Rq. While the ballistic regime SR 2 cx t 2 for 
t <C to, first predicted by Batchelor [33], has been observed 
numerically |34j and experimentally |21j . the existence of 
the Richardson regime (i5R 2 cx t 3 ) is not well established 
(see [35] for a recent review). 

We measured this relative dispersion (SR ■ 5R)(t;R ), 
conditioned on the initial separation i?o, for the three dif- 
ferent types of particles. The statistics of inertial particles 
were expected to be different from that of fluid particles 
since they do not follow the fluid [18]. Figure [3] shows the 
measured data for a particular initial separation normal- 
ized by the Batchelor predictions for tracers. In agreement 
with [3TJ[33], the fluid particles followed almost perfectly 
the Batchclor's regime for times smaller than ~ 0.07to- 
For inertial particles, the t 2 law still held for times smaller 
than ~ 0.07to, but with the prefactor increasing with par- 
ticle density. Different initial separations (for constant 
particle density) showed also Batchelor like scaling, but 
with different prefactors, as discussed later. The increase 
in the prefactor with density reflects that heavier particles 
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Fig. 3: Temporal evolution of the relative separation between 
two particles, for three different types of particle, compared 
with the Batchelor's law. The initial separation Rq is IOO77 
(within the inertial range). 

separated initially faster than fluid particles. 

This last observation, coming from the Lagrangian point 
of view, is inherently connected to the second order Eule- 
rian structure function. The Batchelor regime was de- 
rived from a Taylor expansion of (<5R • <5R)(£;i?o) ~ 
(S r v(Ro) ■ 6 r v(R Q ))t 2 = D 2 (r = R Q )t 2 . As shown in 
Figure HI we measured the Eulerian quantity D 2 (r) in- 
dependently from the particle pair dispersion, which is 
a Lagrangian measurement. The second order velocity 
structure function increased systematically when increas- 
ing the particle density. This is in agreement with the 
observation that heavier particles separated faster in the 
ballistic regime. The structure function D2(r) also in- 
creased with length scale r in the inertial range. There- 
fore the scaling observed and expected for fluid particles 
D% oc (r/77) 2 / 3 applies only approximatively to inertial 
particles. The inset of Fig. 2] shows the evolution of the 
logarithmic slope ^jr 2 ^ °f ^(r) as a function of the 
scale r. The average values of this logarithmic slope, £2, 
for 50 < r/r) < 150 are : & (tracer) = 0.71 ± 0.02 (slightly 
higher than 2/3, possibly due to small intermittency cor- 
rections), 6 (glass) = 0.76±0.02, and £ 2 (steel) = 0.8±0.1 
where the errorbars correspond to the standard deviations 
over the same range of scales. 

To quantify the effect of inertia on relative dispersion of 
particle pairs with initial separations in the inertial range, 
we define the following ratio: 



where X(t) ' corresponds to a time average of X(t) 
for t e [£1; £2]- This quantity T could be viewed as a 
coefficient appearing in the Batchelor regime for inertial 
particles: (SR-6R) = ^-T (R a ) R^) 2 . Note that T > 1 



Fig. 4: Second order Eulerian structure function D2(r) = 
(5 r v ■ <5 r v) normalized by 11 ^' 2 (er) 2,/3 the inertial range value 
of D2{r) expected for fluid particles. The inset represent the 
logarithmic slope of D2 (r) for the three different types of par- 
ticles, the horizontal dashed lines represent the average of this 
quantities for 50 < rjr\ < 150. The symbols are the same as 
the one used in Fig. [3] 

means that the separation velocities between particles are 
larger than that for fluid tracers. This ratio, obtained 
separately from the Lagrangian measurement of relative 
dispersion and from Eulerian measurement of the second 
order velocity structure functions, is shown in Figure [5] 
In the inertial range (Rq j§ 50?/) T is above one and 
increases with scale, which confirms that in this range 
of scales, the inertial particles separated faster than the 
fluid particles. Moreover, the independent Eulerian and 
Lagrangian measurements of T coincide almost perfectly. 
In the experimentally accessible range of inertial scales 
50 < Ro/r] < 150, as shown in Fig. [5j T is mostly greater 
than unity, and is increasing with Rq. This reflects that 
heavy particles separate faster than tracer particles and 
that the scaling of the second order velocity structure func- 
tion is affected by the particle inertia. In their numerical 
simulations, Salazar & Collins [3Zj found a similar effect 
even if their inertial range is much smaller because of a 
small Reynolds number (R\ = 120). 

At sufficiently large separations (on the order of the in- 
tegral scale L ~ 10 3 ?/), the velocity spatial correlation 
vanishes, therefore T tends to the ratio of the average 
particle kinetic energy per unit mass k p = |(v • v) to 
that of the fluid kf = |(u • u). In our measurements, 
surprisingly, this ratio is bigger than unity and varies 
non-monotonically with increasing particle density p p : 
kp(p p = i)/kf = 1.063 and k p (p p = 7.8)/fc/ = 1.05. In a 
recent numerical work, Salazar & Collins [37| observed the 
same effect using a point-particle model: With increasing 
St, this ratio first increased to St J$ 0.25 and then fell bel- 
low unity. Moreover, by computing the fluid kinetic energy 
along the inertial particles trajectories they demonstrated 
clearly, that for St < 0.5 the non-monotonic behaviour of 
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Fig. 5: Ratio T(Ro) (see text). The colors correspond to the 
two types of inertial particles considered, red for the glass parti- 
cles (pp/pf — 4) and black for the steel particles (pp/pf — 7.8). 
The solid-lines are the Eulerian measurements. Symbols with 
errorbars are the Lagrangian measurements. The dashed-lines 
correspond to the model developed in the text. 

k p with St is due to the uneven sampling of the flow by 
inertial particles. 

As an attempt to understand analytically the observed 
increase of relative velocity with particle inertia, we can 
start with the simplest equation of motion for the inertial 
particles: 

e?v „Du 1 . „. 

— =/3— + - u-v)+ l-/3 g, 7 

at Dt t p 

where v and u are the velocities of the particles and the 
fluid at the particle position, the modified density ratio 
(3 = 3pf/(2p p + Pf) takes into account the added mass 
effect, t p = dp/ 12 f3 v is the particle viscous relaxation time, 
g is the gravity vector and the operators d/dt and D/Dt 
indicate derivatives following an inertial particle or a fluid 
element, respectively. Eq. ([7]) is the simplest form of the 
equation derived for inertial particles (see, e.g. [THIUS]) 
and has been used before to study the effect of particle 
inertia [38] . Using this equation, one can solve for particle 
velocity in the limit of small r p (i.e., St -C 1): 

v=u + r p (l-/3)(g-a) + (r p 2 ) (8) 

With this we evaluate the Eulerian structure function of 
the inertial particles in terms of statistical quantities of 
the flow field: 

(5 r v ■ 5 r v)= (S r u ■ 5 r u) + 2(/3 — l)r p {<5 r u • S r a) + o(Tp), 

(9) 

where a = Du/Dt is the fluid acceleration. We can com- 
pare directly the results of this model to the measurements 
by computing T(r) (see equation : 

(<Vu-o r u) IIC2 \r] J 



This expression is only valid for length scales in the inertial 
range (r ^ 50?y) since we used the relation (6 r u-6 r a) = —2e 
as given in Eq. (g]). For heavy particles (p p > pf), the 
density ratio /3 < 1. Our simple equation predicts that 
r > 1, i.e., heavier particles possess larger relative ve- 
locities and hence separate faster. This is in qualitative 
agreement with the experimental observation. However, 
as shown in Figure the quantitative agreement between 
the simple model and the measurements is not satisfac- 
tory. In particular, the simple model shows that the rel- 
ative increase of relative velocity between inertial parti- 
cles is proportional to a "scale-dependent Stokes num- 
ber" St r = St(rj/r) 2 / 3 = Tp/r r , which compares the par- 
ticle time scale to the life time of an eddy of scale r, 
t> = (r 2 /e) 1 / 3 and hence decreases with scale r, while 
the measurements shows that effect increases with scale 
in the range accessible experimentally. 

We now discuss the reasons why this simple model 
does not agree quantitatively with the measurements. At 
first sight, it might come by neglecting terms of order t 2 
and higher. However, as one can see from Figure [5] if 
this was the case, it would mean that the higher order 
terms are actually the leading terms, which is unlikely. 
The second possible explanation comes from the fact that 
in the last step of Eq. (fTU)). we implicitly assumed that 
inertial particles distribute evenly in space. This assump- 
tion ignored the so called "sampling effect" of inertial 
particles, as it has been shown that heavy particles are 
ejected by intense vortices and tend to accumulate in 
low-vorticity/high-strain regions [H1I32]- Therefore, these 
heavy particles might experience a higher "dissipation 
rate". As a first order correction to this effect, we 
replaced the term (S r u ■ S r a) in Eq. © by the measured 
value of (<5 r v • <5 r a p ), where a p was the acceleration of 
inertial particles. The calculated T(r) still does not agree 
with measurements, especially its dependence on scale 
r. Using numerical simulations, Salazar & Collins |37| 
showed that Eq. [5] led to T(r) that agrees reasonably well 
with their DNS if the quantities on the right hand side 
of Eq. [HI were conditioned on the particle trajectories. In 
order to further identify, whether sampling or filtering 
|37j is responsible for the apparent discrepancy between 
the model and the data, we estimated with Eq. [5] the 
fluid velocity at the particle position , up to order r 2 
(using a p = a + o(t p )). For all three Stokes numbers the 
second order velocity structure functions of this estimated 
fluid velocity field showed the same logarithmic slopes as 
those of the measured particle velocity field. Thus, the 
lack of agreement between the experimental data and the 
simple model, Eq. I1Q[ most likely roots in the uneven 
sampling of the flow field by inertial particles and not 
in the filtering of the velocity field by inertial particles. 
At the Stokes numbers observed the non-differentiablc 
particle velocity field (caustics) can also contribute as 
discussed by Salazar & Collins [57], however, we believe 
this to be a less important effect. 
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We summarize our experimental results as follows. For 
heavy, Kolmogorov-sized particles we have observed a 
ballistic, or Batchelor regime for times t ;$ O.Uq with 
to determined by the turbulence energy dissipation rate 
and the initial separation between particle pairs. In this 
regime, these heavy, inertial particles separated faster 
than fluid tracers, i.e., the relative velocity between parti- 
cles increased with particle inertia. This was also reflected 
by the measured Eulerian velocity structure functions. 
In the inertial range, the logarithmic slope of the second 
order velocity structure function increased with the parti- 
cle density. We also observed a non-monotonic variation 
of the RMS velocity with Stokes number, which, in this 
range of Stokes numbers, Salazar & Collins [37] attributed 
to the preferential sampling of the inertial particles. The 
observed effect of inertia on particle dynamics could not 
be captured by a model based on the simplified equation 
of motion for idealized point-like particles. The inaccu- 
racy of the model may originate from the simplified model 
assumptions that do not capture the preferential sampling 
of inertial particles |40j . In addition, a point-particle 
model might not be sufficient to describe the dynamics 
of the Kolmogorv size particles studied here. Finally we 
note that the increase of relative velocity between inertial 
particles with their density will result in an increase 
of the collision rate between inertial particles, which 
has important consequences in problems such as rain 
formation in warm clouds. To understand these issues 
quantitatively, we would need to spatially resolve the flow 
at the particle scale. This is currently under investigation. 
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